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Abstract. Considered in this paper the Maxwell-Bloch (MB) equations became known after 
Lamb U-0]. In [5] Ablowitz, Kaup and Newell proposed the inverse scattering transform (1ST) 
to the Maxwell-Bloch equations for studying a physical phenomenon known as the self-induced 
transparency. A description of general solutions to the MB equations and their classification 
was done in [6] by Gabitov, Zakharov and Mikhailov. In particular, they gave an approximate 
solution of the mixed problem to the MB equations in the domain x,t £ (0, L) X (0, oo) and, on 
this bases, a description of the phenomenon of superfiuorescence. It was emphasized in [6] that 
the 1ST method is non-adopted for the mixed problem. Authors of the mentioned papers have 
developed the 1ST method in the form of the Marchenko integral equations. We propose another 
approach for solving the mixed problem to the MB equations in the quarter plane. We use 
a simultaneous spectral analysis of the both Lax operators and matrix Riemann-Hilbert (RH) 
problems. First, we introduce appropriate compatible solutions of the corresponding Ablowitz- 
Kaup-Newell-Segur (AKNS) equations and then we suggest such a matrix RH problem which 
corresponds to the mixed problem for MB equations. Second, we generalize this matrix RH 
problem, prove a unique solvability of the new RH problem and show that the RH problem (after 
a specialization of jump matrix) generates the MB equations. As a result we obtain solutions 
defined on the whole line and studied in [5] and [6], solutions to the mixed problem studied 
below in this paper and solutions with a periodic (finite-gap) boundary conditions. The kind 
of solution is defined by the specialization of conjugation contour and jump matrix. Suggested 
matrix RH problems will be useful for studying the long time/long distance (x £ K+) asymptotic 
behavior of solutions to the MB equations using the Deift-Zhou method of steepest decent. 



1. Introduction 

The Maxwell-Bloch equations arise in different physical problems. Most significant applications 
of this system deal with the problem of the propagation of an electromagnetic wave (ultrashort 
optical pulse) in a resonant medium with distributed two- level atoms. In particular, there are 
the problem of self-induced transparency, the laser problems of quantum amplifier and super- 
fluorescence. The system of the Maxwell-Bloch (MB) equations can be written in the form 

(1.1) £ t + £ x =(p), 

(1.2) Pt + 2i\p =M£, 

(1.3) N t =- l -{£*p + £p*). 

Here £ = £(t,x) is the complex electric field envelope, so that the field in the resonant medium is 

V(t,x) =£{t,x)e in ( x - t] +£*(t,x)e- in{x - t \ 
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Subindexes mean partial derivatives in t and x, and * means a complex conjugation. M = Af(t, x, A) 
and p = p(t, x, A) are entries of the density matrix of a quantum two-level atom subsystem. The 
parameter A is the deviation of the transition frequency of the given two-level atom from the mean 
frequency £1. The angular brackets mean averaging 

oo 

(p) = n J n(X)p(t,x,X)d\ 

— oo 

with the given "weight" function n(A), such that 

oo 

J n(A)dA = ±l. 

— oo 

The weight function n(A) characterizes the inhomogeneous broadening, i.e. a form of the spectral 
line. It is the difference between the initial populations of the upper and lower levels. If n(X) > 0, 
then an unstable medium is considered (the so-called quantum laser amplifier). If n(A) < 0, then 
a stable medium is considered (the so-called attenuator). 

The mixed problem to the MB equations is defined by initial and boundary conditions 

(1.4) £(0, x) =£q(x), < x < L < oo 

(1.5) p(0, x, A) —po(x, A), 

(1.6) Af(0,x,\) =Af (x,\), 

(1.7) £(t,0) =£ in (t), 0<t<oo. 

We assume that entering in the medium < x < L pulse £i n {t) is smooth and fast decreasing. 
Initial functions £q(x) and po{x, A) are also smooth for < x < L and A e M. Function Nq{x, A) 
is defined by po(x, A): 



(1.8) M)(i,A) = VH po(x,X)\ 2 , 0<x<L<oo. 

Here we chose the positive branch of the square root. To find (|1.8p we use the second and the third 
equations of (|l.l[) - (|1.3p . They give 

^- t (N 2 (t,x,X) + \p(t,x,X)\ 2 ) = 0, 

and one can put 

Af 2 (t,x,X) + |p(t,a;,A)| 2 = 1. 
In what follows we also put 57 = 1 in (|1.4[) . 

The Lax pair for the Maxwell-Bloch system was first found in [5] by using results of [1]- [4] (see 
also [6] and [7]). It was shown that (|l.l[) - (|1.3j) are equivalent to the over determined linear system, 
known as Ablowitz-Kaup-Newell-Segur (AKNS) equations: 

(1.9) Wt + iXcr^w = — H(t t x)w, 

(1.10) w x — 1X03W + iG(t, x, X)w = H(t, x)w, 
where 173, H (t, x), G(t, x, A) are the matrices 

(I \ 1 / £{t,x) 

OO 
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and the symbol p. v. denotes the Cauchy principal value integral. Differential equations ()1.9j) and 
()1.10P are compatible if and only if £ (t, x), p(t, x, A) and Af(t, x, A) satisfy equations (|l.ljl - (|1.3jl (cf. 
[7]). As shown in [3], p(t,x,X) and Af(t,x,X) are related to the fundamental solutions of (|1.9[) . 
Indeed, let <P(t, x, X) be a solution of (jl.9[> such that det <P(t, x,X) = 1. Then it is easy to show that 
F(t, x, X) = <P(t, x, X)a3<I>i(t, x, A), where & is the Hermitian-conjugated to <P, satisfies equation 

(1.11) F t + [iXa 3 + H,F]= , F{t,x,X)=(£f>Z>X> ^ x %) 

which is a matrix representation of equations (|1.2j) - (|1.3j) . It is well-known that the fundamental 
solution of equation (ll.9|) is entire in A. Therefore the matrix F(t, x, A) = <P(t, x, X)a^<P^(t, x, A) is 
smooth in A € R if initial matrix F(0, x, A) is smooth too. 

In what follows we will use the x ± -equations (upper/lower bank equations instead of x-equation): 

(1.12) w x — iXa^w + iG±(t, x, X)w = H(t, x)w, 
where 

oo oo 

„ , \ 1 f F(t,x,s)n(s) , 1 f F(t,x,s)n(s) , , 7ri , , . ,,. 

G±(t,x,X) = - J )_ X ^> Q J ds=p.v.- J K s _' x y > d S ±-F(t,x,X)n(X). 

— oo — oo 

Thus we have new Lax pairs (t- and x + -equations and t- and x~ -equations) for the MB equations. 
Equations (|1.9[) and (|1 . 12|) (as well as (| 1 .9[) and (|1.10|l ) are compatible if and only if £ (t, x), p(t, x, A) 
and Af(t, x, A) satisfy equations (|1.1|) - (|1.3|) . This known fact is proved below (see section 2). 

Main goal of the paper is to show that the inverse scattering transform is applicable to the mixed 
problem for the MB equations using a simultaneous spectral analysis of the compatible differential 
equations (11.91) and (|1.12p . We develop the 1ST method in the form of matrix Riemann-Hilbcrt 
problem in the complex plane and give an integral representation for £ (t, x) through the solution 
of a singular integral equation, which is equivalent to the matrix RH problem. To produce this 
RH problem we introduce appropriate compatible solutions of the corresponding AKNS equations 
and input spectral (scattering) data defined through given initial and boundary conditions for MB 
equations. Then we suggest such a matrix RH problem which corresponds to the mixed problem for 
MB equations. Further, we give a formulation of more general matrix RH problem, prove a unique 
solvability of the new RH problem. After a specialization of jump matrix we prove that the RH 
problem generates the AKNS linear equations for the MB equations with a given inhomogeneous 
broadening. Thus this RH problem generates different solutions to the MB equations. There are 
solutions defined on the whole line and studied in [5] and [5], solutions to the mixed problem (|1.4I) - 
(|1.7p with vanishing boundary conditions and periodic (finite-gap) boundary conditions, etc. The 
kind of solution is defined by the specialization of conjugation contour and jump matrix. Suggested 
matrix RH problems will be useful for studying the long time/long distance (x € R + ) asymptotic 
behavior of solutions to the MB equations using the nonlinear method of steepest decent. 

Our approach differs from that was proposed in [TU] for the Goursat problem to the MB equa- 
tions. We develop the method initially proposed in [TT]-[T3] and also in [15]- [IB]- This method uses 
the simultaneous spectral analysis of the Lax operators. We formulate such a matrix RH problem 
which allow to prove that the mixed problem is linearizable completely. 

We will use the following notations: if A denotes matrix A = (A[l], A[2]), then vectors A[l], 
A[2] denote the first and the second columns of matrix A. We also set [A, B] = AB — BA. 
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2. Compatible solutions of the AKNS equations 

Let 

(2.1) W t = U(t,x,X)W, 

(2.2) W x = V(t,x,X)W 

be a system of matrix differential equations. Equations (|2.1j) . (|2.2p are compatible if W tx = W xt 
for any solution of these equations, i.e. the following condition holds : 

(2.3) U x - V t + UV - VU = 0. 

Lemma 2.1. Let U(t,x,X) and V(t,x, A) are defined and smooth for all t, x, A € R and satisfy 
condition (|2.3p . Let W^tjX, A) satisfy t-equation (|2.1|) /or aZ/ x, and Zet Vy(to,x, A) satisfy x- 
equation (|2.2| /or some t = to, ^0 < °°- TTien W^(t, x, A) satisfies x-equation for all t. The same 
result is true, if one changes t into x and vice versa. 

Proof. Let W(t, x, A) = W x — V(t, x, X)W . Then W(t, x, A) solves equation flU]). Indeed 

W t = U(t, x, X)W + (U x -V t + [U, V})W = U(t, x, X)W. 

Matrices W , W are solutions of the same equation. Therefore they are linear dependent, i.e. 
W{t, x, A) = W(x, t, X)C(x, A). Since W(t ,x, A) = 0, then C(x. A) and VK(t, x, A) are identically 
equal to zero. Lemma |2. H is proved. □ 

Due to (fl~9l) . (fETO) and (fl~T2|) we have three AKNS systems. The t-equation ([H)]) is defined by 
matrix U (t, x, A) = —iXcr^ — H(t, x) while the x-equation (|1.10|) and upper/lower bank x ± -equations 
(|1.12|) are defined by matrices: 

V(t, x, A) =iAcr 3 + H{t, x) - iG(t, x, A), 
V ± (t,x,X) =V(t,x,X)±'^^-F{t,x,X). 

The compatibility condition of t-equation (|1.9I) and x-equation f| 1 . 10[) takes the form of (|2.3[) . The 
compatibility condition of the t-equation (jl.9[) and x ± -equations (|1.12|) is as follows: 

[4 - Vt + [U, V±] = U X -V t + [U, V] T ^^{F t + [iAaa + H, F}) = 0. 

If F t + [iAo-3 + H , F] = (in view of eq. (jl.lll) ) , then the new compatibility conditions coincides with 
Conversely, if U x -V t +[U, V] = 0, then F t +[iXa 3 +H, F] = as well as U x -Vt + [U, V*] = 0. 
Each of them gives the same MB equations. The upper/lower bank equation (|1.12[) (in a contrast 
with x-equation (jl.lpp ) will allow to obtain such vector-solutions of the AKNS equations, which 
have analytic continuation to the upper and lower complex z-planc (z = A + izy). We will use them 
for a construction of the matrix Riemann-Hilbert problem, which gives a solution of the mixed 
problem to the Maxwell-Bloch equations. 

Theorem 2.1. Let £(t,x), Af(t,x,X), p(t,x, A) be a smooth solution to the MB equations U.l\) - 
11. Sty defined for all t,x € R+ and X € K. Let £(t, 0) is fast decreasing as t — > +oo and let 
p(0, x, A) = for x > L or is also fast decreasing as x — > +oo ( if L = oo). Then there exists two 
pairs of compatible solutions Y^ft, x, A) and ^ ± (t, x, A) of the AKNS equations U.9\) . \1.12}) such 
that 

(2.4) Y ± (t,x,X) = W ± (t,x,X)$(t,X), 



(2.5) 



Z ± (t, x, A) = *(t, x, A)w ± (x, A). 
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Here W^(t, x, A) satisfies x^ -equation S1.12\) for all t, W(t, 0, A) = /, and <E>(t, A) satisfies t- 
equation M.9\) with x — under the Jost initial condition: 

lim $(t, A)e iAtCT3 = /. 

Function ^(t,x,A) satisfies t-equation $1 . 9\) under initial condition ^(0,x,A) = J ; and zt^l^A) 
satisfies x ± -equation \1.12\) with t = under initial condition w (L, A) = e iiL77± ^°" 3 ? where 



oo 

(2.6) »?±(A) = A-i | 



n(s)ds 
s - A i iO ' 



— oo 



If L = oo then the condition at x = L changes with the Jost initial condition: 

lim w ± {x 1 X)e-' 171 ^ xa:t = I. 

x— >oo 

Proof. If matrices W^i, a;, A), A) and ^(i, a;, A), w ± (x 7 A) do exist then, due to Lemma 2.1, the 
products M /± (i, x, A)<J>(t, A) and \I/(t, a;, A) w 1 * 1 (x, A) are compatible solutions of the AKNS equations 
(|1.9|) . (|1.12l) . The existence of these matrices will be given in the next Lemmas. □ 



*(U) 



(t,x) 



W ^ X ' A )^(t,X,A) 



w±(x,A) L X 



Lemma 2.2. Let £(t,0) = £\{t) be smooth and fast decreasing as t — > oo. Then the Jost solution 
<&(t, A) has an integral representation: 



(2.7) A) = e _lAt<T3 + y Af(i, r)e~ lAr<T3 dr, Im A = 0. 

The kernel K(t,r) satisfies the symmetry condition K*(t,r) = AK(t, r)A _1 with matrix A 
and [0*3, K(t,t)] = o~3H(t,0). The kernel K(t,r) is smooth and fast decreasing as t + r 



-1 
oo. 

The proof of this Lemma is well-known (cf. [19). Due to this integral representation, vector- 
columns <&[l](i, A) and <&[2] (t, A) of the matrix $(i, A) = ($[l](i, A), <I>[2](^ A)) have analytic con- 
tinuations $[l](i, z) and $[2](i, z) to the lower and upper half planes of the complex z-planc 
respectively. The next asymptotic formulae are valid: 

• $[!](*, z)e izt = ) +0{z~ 1 ), lmz<0, z ->• oo, 



• $[2](t, z)e~ izt = n +0(z- 1 ), lmz>0 z ->• oo. 
They are a simple consequence of integral representation (|2.7p . 
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Lemma 2.3. Let n(X) be a Holder function for A £ 1 and let £(t,x), Af(t,x,X), p{t,x,X) be 
smooth functions for t,A £ R and x G R+. Then the solution W (t, x, A) can be represented in the 
form: 

(2.8) W ± (t,x,\)=e iXx ^ X ± (t,x,\), 
where ^^(t,^, A) is the unique solution of the Volterra integral equation 

X 

(2.9) X ± {t,x,X)=I + J e- iX y^(H(t 1 y)-iG ± (t,y,X))e iX y^ X ± (t,y,X)dy, 

n 

The solutions W^(t,x, X) are smooth in t and x, related each other by the formula: 



W-(t,x,X)=a 2 W + *(t,x,X)a 2 , a 2 



i 

-i 



and have an analytic continuations W(t,x,z) to the upper C+ and W(t, x, z) to the lower C_ 
complex z-plane respectively. Moreover, W(t, x, z)e~ lzx is continuous, bounded in C_ UR and have 
the following asymptotics: 



W{t,x,z)e-^=(\ e _L) +0(z- 1 ), 



Im z < 0, z — > oo, 



and matrix W(t,x, z)e lzx is continuous, bounded in C+ Ul and have the following asymptotics: 

W(t,x,z)e izx = (^ 2l * X °^j+0(z- 1 ), lmz>0. z -> oo. 

Proof. The solvability of the Volterra integral equation (|2.9p and the smoothness of the solution 
with respect to t and x can be easily proved by using the method of successive approximations. So- 
lutions W (t, x, X) are not independent. One can easily verify that W~(t, x, X) = o-2W + *(t, x, X)a 2 . 
Since G±(t, x, A) has an analytic continuation to the domain C±, the matrix W (t,x, X) has also 
an analytic continuation for z = X + \v G C± which we denote as W(t, x, z) for z G C+ and 
W(t,x,z) for z G C_. Thus the matrix W(t, x, z)e~ lzx (W(t, x, z)e lzx ) is analytic, continuous and 
bounded in C_ U R (C+ U R) and has the above asymptotic behavior. □ 



The previous lemmas and formulas (|2.4|) . (|2.7|) . (|2.8|) . (|2.9[) imply the following properties of the 
matrix Y ± (t,x,X) = {Y ± [l]{t, x, X) Y ± [2](t, x, A)): 

(1) Y (t,x, A) satisfy the t— and a; 1 * 1 — equations; 

(2) dety ± (f,a;,A) = 1, Ael; 

(3) the map (x, t) i — > Y (t, x, A) is smooth in t and x\ 

(4) vector column Y~ [l](t, x, A) has analytic continuation x, z) for z G C_, which is 
continuous and bounded in z G C_ U R and 

Y[l]{t,x,z)e izt -' lzx = Q +0(z~ 1 ), z^oo; 

(5) vector column Y~ + [2](i, a;, A) has analytic continuation y[2](<, z) for z G C+, which is 
continuous and bounded in z G C+ U M and 



r[2](t,x,z)e 



-\zt-\-'\zx I ' ' 



j)+0(z-!), 
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Lemma 2.4. Let £(t,x) be smooth functions in the domain of its definition. Then the function 
ty(t,x,\) has an integral representation: 

(2.10) V(t,x,\) = e~ ixta3 + f L(t,T,x)e~ iXT,73 dT. 



The kernel L(t,T,x) is smooth, satisfies the symmetry condition L*(t,T,x) = AL(t,T, x)A 1 with 
1> 
-1 0, 



matrix A = ( 3.1 and [173, L(t,t, x)] = H(t, x)as 



The proof is well-known (sf.[15|). This lemma gives that 9(t, x, A) has an analytic continuation 
ty(t, x, z) for all z £ C. Moreover ty(t, x, z)e lzt is continuous, bounded in C+ U K and 

(2.11) 9(t,x,z)e^=(l e L)+0(0. l«|->«>, 
and z)e _lzt is continuous, bounded in C_ U K and 

(2.12) *(t s ar,2!)e- fa *= ^"q"* ^+0(0. 1*1 "> °°- 

If H(t,x) = and F(t,x,X) = 03, then adequation has the exact solution e 13 " 7 *'*^ 3 , where 
J?±(A) is defined in (12. 6|) . These functions are the boundary values of the sectionally analytic 
function 



(2.13) r>{z)=z-\ f ^ds 



—00 



for z G C±. Across the real line it has the jump: 



r 7 (A + iO)-7 ? (A-iO) = -yn(A) Ae 



Furthermore, for z = A 



If n(A) < then sign Im n(z) = signlmz. In the case n(X) > 0, 

Im?7(z) 



Im z 
Imr](z) 



>0 if I(A,i/)<1; 
<0 if I(A,i/)>1. 



Im z 

For n(A) = 5(A) this integral equals to 1/4(A 2 + v 2 ) = l/|2z| 2 and equation /(A, v) = 1 defines the 
circle |z| = 1/2. In a general position when n(A) > 0, there exists a curve, which is the boundary 
of a domain D containing the origin of the complex z-plane. This curve is defined by equation: 
Imrj(z) — when Imz ^ 0. It is symmetric with respect to real A-axis, because 7]*(z*) — r](z). We 
denote it as 7 U 7, where 7 lies in C+, and 7 in C_. Thus there is such a domain D = D + U D~ 
that 

-signlmz, z E D = D + U D~ , 
signlmz, z € (C+ \ D+) U (C_ \ D~). 



signlm77(z) 
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Domain D may be bounded and unbounded as well. Indeed, let 

n(A) 



0, |A|>£. 

Then the curve 7 (Im77(z) = 0,Im2: 7^ 0) is described by equation: 



v ' V tan8e^ 
where the square root is positive. If e is positive and sufficiently small then there exists 5 > such 
that 



4 + 5 w V 4' 
i.e. finite curve 7 together with interval [A_,A + ] (A± = i-^/1/4 + e 2 ) bound domain D + . An 
example of the unbounded domain D gives a physical and simple model with the Lorentzian line 
shape: 

Then 

1 



r,{z) = z-- —dX 



-■4(* + i0' ZGC+ 
Z+ 4(^i0' * eC " 



and curve 7 is defined by equation: 



This equation yields that A±(0) = ±00, i.e. (A_,A + ) = M, and z/ max = (y/l + 1 2 — l)/2. Thus 
domain D + is unbounded along the real line because \{v) is unbounded when v — > 0. Domain D~ 
being complex conjugated to D + is unbounded also. 
If -F(0, £, A) = 1T3 for x > L we put 

JS^foA) =i7(0,x)-i(G±(0,x,A)- 5± (A)cr 3 ), < x < L, 



fl ± (A) = \ J 



n(s)ds 
s - A =F i0 ' 



Lemma 2.5. Let n(A) 6e a Holder function and fast decreasing as A — > ±00. Lei £(0,x) = £q{x), 
Af(0, x, A) = Afo(x, X), p(0, x, A) = /Oo( x i A) be smooth for AeK, 0<i<I and fast decreasing as 
x — > 00 w/ien L = 00. Then w ± (x, A) can 6e represented in the form: 

(2.14) ur^s.A) - w ± (x,A)e i " ± ( A ) 2:<T3 , 

w/iere w ± (x, A) is i/ie unique solution of the Volterra integral equation 

L 

(2.15) w ± (x, A) = J - J ^ ± ^ x -^'"H ± (y, A)w ± (y, AK^^-^dy. 

The solutions w ± (x, X) are smooth in t and x, related each other by the formula: 

w~(x, A) = a 2 w + *(x, A)o2- 
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If n(X) < then the first column u> + [l](x, A) has an analytic continuation w[l](x,z) to the 
upper C+ complex half-plane and the second column w~[2](x, A) has an analytic continuation 
w[2](x, z) to the lower C_ half-plane. The corresponding analytic vector-column w[l](x, z)e^ lr, ^ x 
and W[2](x, z)e lrt ^ x are continuous in C± UR and have the following asymptotics: 



w[l]{x, z)e-'"'^ x = f ^ + Oiz- 1 ), Imr)(z) > 0, 



'2](x,z)e iT >W x = ("] +O0" 1 ), Imr?O)<0, z 




1 

If n(X) > i/ien t/ie first column w + [l](x, A) and £/ie second column w~[2](x,X) /iaue analytic 
continuations into domains C + \ D + and C_ \ D~ respectively with the same asymptotic behav- 
ior, while the first column w _ [l](x, A) and the second column w + [2](x, A) have analytic (bounded) 
continuations w[l]{x , z) and w[2](x, z) to the domains D~ and D + respectively. 

The proof can be done by the same way as for Lemma 2.3. 

Lemmas 2.4, 2.5 and formulas ([23)1. (|2~TU|) . ([2~TTT) . ([2~T2l . (|2Tl"3)) . (|2~T51) imply the following 
properties of the matrix Z ± (t, x, A) = (Z ± [l](i, x, A), Z ± [2](i, i, A)): 

(1) Z (t,x, A) satisfy the i— and x 1 * 1 — equations; 

(2) det ^(i, x, A) = 1, Ael; 

(3) the map (x, t) i — > Z ± (t, x, A) is smooth in £ and x; 

(4) vector column x, A) has an analytic continuation Z[l](f , x, z) for z € C + if n(X) < 
(z e C+ \ D + if n(A) > 0), which is continuous in the closure of C+ (C+ \ D + ) and 



Z[l](t, x, z)e 



izt — irj(z)x 



0{z 



(5) vector column Z [2] (i, x, A) has an analytic continuation Z[2] (i, x, z) for z e C_ if n(X) < 
(z e C- \ D~ if n(A) > 0), which is continuous in the closure of C_ (C_ \ D~) and 



Z[2](t,x,z)e 



-izt-\-ir](z)x 



(J)+o(^ 1 ), 



If n(A) > then 

(6) vector column Z~[l](t,x, A) has an analytic continuation Z[l](t, x, z) to the domain D~, 
and vector column Z + [2](t,x, A) has an analytic continuation Z[2](t, x, z) to the domain 
D + , where the both analytic continuations are continuous up to the boundary of D + and 
D~ respectively. 

To prove (4) and (5) we use asymptotic relations (12. lip . (|2.12[) and following one (where signs ± 
are omitted for a convenience) 



w(x, A) 



e ixri(\) 






e -ixr)(X) 



X u(x, X)e ixr >W X i2(x,X)e- ix ^ 
X 2i(x,X)e ixv(X) X 22(x,X)e^ x,l{X) 



Here X ii(x,A), X2i{ x , A) = 0(A~ 1 )+0(A~ 1 e 2iiA ) and Xi2(x, A), X22(x, A) = 0(X~ 1 )+0(X~ 1 e~ 2iLX ) 
asA^ ±oo (0 < x < L < oo). If L 
first vector column, we have 



-l e 2iLA>j 

oo then Xij(x, A) = 0(A _1 ) as A 



-2iLA-j 

±oo. Indeed, for the 



Z[\\ (t, x, z )e^-^> ~((\ +0(z- 



For the second vector column there is 

—2izt q 

o 1 



Z[2](t,x,z)e 



-\zt- s r'\r){z)x 



+ 0(z~ 1 ) 



X2l(x,z) 

Xl2(x,z) 
X22(X,Z) 



0(z" 



+ 0(z" 1 ). 
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Due to analytic continuations of the first vector column to the upper half plane and second 
vector column to the lower half plane we have (xu(i,z), X2i(x,z)) = f° r Imz > and 

(X2i(x, z), X22(x, z)) = 0(z~ 1 ) for Imz < as z — >• oo. 

Pairs of vectors (Y[l](f, x, z), Z[2](t, x, z)) and (Z[l](t,x, z), Y[2](t,x, z)) will alow to define 
below suitable matrix Riemann-Hilbert problem in the case n(X) < 0. In the case n(A) > we will 
use additionally more vectors (Z[l](t,x,z), Y[l](t,x,z)) and (Y[2](t,x, z), Z[2](t,x,z)). 

At the end of this section we consider the problem on whole i-line studied in [5] and [B] , where 
the Marchenko integral equations were used. We would like to formulate below (in section 4) the 
corresponding matrix RH problem. Let £ (t, x) is vanishing as t — > — oo in such a way that there 
exists the Jost solution 

t 

#(t,x, A) = e- iXta3 + / L(t,T,x)e' iXTCT3 dT 



of i-equation (jl.9p . Let Af(t,x,X) — 03 and p(t, x, A) are also vanishing as t — > —00. In this case 
we use the following pair of compatible solutions: Y (i 5 x, A) are the same as above, while 

Z ± (t,x,X) = ^(t,x,\)e ixn±Wa3 . 

Due to the Lemma 2.1 Z ± (t, x, A) is compatible solution of (|1.9j) and (|1.12j) because, when 
Af(t,x,\) — 03 and p(t,x,X) — 0, x ± -equation takes the form E x — i?7± ( A)<r 3 i? and we put 
E ± (x, A) = e lxr,± ^' T3 . Further, vector column Z + [l](t, x, A) has an analytic continuation Z[l)(t, x, z) 
for z £ C+, and vector column Z~[2](t, x, A) has an analytic continuation Z[2](t,x, z) for z G C_. 
They satisfy the properties (4) and (5) above for any choice of the sign of the function n(X). 

3. Transition matrices and spectral functions 

Now we consider a dependence between solutions of the AKNS equations, corresponding tran- 
sition matrices and spectral functions. First of all, since Y ± (t,x,X) and Z ± (t, x, A) are solutions 
of the t— and a; 1 * 1 — equations (|1.9p . (|1.12p . they are linear dependent. So there exist transition 
matrices T ± (A), independent of x and t, such that 

(3.1) Y ± (t,x,X) = Z ± (t,x,X)T ± (X), ImA = 0. 

They are equal to 

T^A) = (w ± (0,A))- 1 $(0,A), ImA = 
and have the following structure: 

T±m /a±(A) 6±(A) 
T (A) = {-t(X) a±(X) 
Matrices w ± (0, A) have the similar structure 

and define the spectral functions of x ± -equation. Since coefficients of a; -equation are taken 
for the fixed time t = 0, the spectral functions are uniquely defined by given initial functions 
£(0,x), p(0,x,X) and 7V(0, x, A). They are not independent because the solutions w ± (x,X) of 
the x ± -equations are related each other by w~(x,X) = <T2W + *(x, X)a2, that yields the following 
reductions: 



a±(A) =a =F *(A), 



r(A)=-^*(A). 
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The matrix 



mX '-{-B(X) A(X)J 



gives the spectral functions of the ^-equation with x = 0. They are uniquely defined by the 
boundary condition £(t, 0). Analogous reduction condition A) = CT2<& + *(£, X)a 2 means 

A(A)=A*(A), B(X)=B*(\). 

Finally we have that the spectral functions, generated by the transition matrix T ± (A), are defined 
by the initial and boundary conditions. This matrix satisfies the analogous reduction: 

T-(A) =<r 2 T+*(X)a 2 , 

which is equivalent to the following ones: 

a±(A) =a T *(A), 

(3.2) b ± (X)=b^*(X). 

Due to A.S.Fokas [TT] we call all of these functions as the spectral functions. 
Entries of the matrix (A) are 

a ± (X)= dct[F ± [1] (t, x, A), Z ± [2] (t, x, A)] 

fe ± (A) = det[F ± [1] (t, x, X),Z ± [1] (t, x, A)] 

6 ± ( A) = det [Y ± [2](t,x,X),Z ± [2](t,x,X)} 

a ± (X) = det^lKt, x, A), Y ± [2](t, x, A)]. 

If n(A) < these relations show that a + (A) has an analytic continuation a(z) for z G C+, and 

a (X) has an analytic continuation a(z) for z 6 C_ (S(z) = a*(z*)). Functions fo ± (A) and 6 (A) 
are defined for A e R only. Determinant of T ± (A) = 1 for ImA = 0. The spectral functions have 
the following asymptotics: 

a(z) =1 + 0(z _1 ), lmz<0, \z\ -> oo. 

6 ± (A) =0(A- 1 ), A^±oo, 

a(z) =1 + 0(z _1 ), Imz > 0, |s| oo. 

If a(z) have zeroes Zj S C+ (j = 1, 2, ....,p) then 

a(zj) =6eb[Z[l](t,x,Zj),Y[2](t,x,Zj)] =0. 

Hence vector-columns of the determinant are linear dependent: 

(3.3) Y[2](t,x,z j )=-y j Z[l](t,x,z j ), 7i = = 

a l z iJ P\ z j) 

At the conjugated points z* 6 C_ (j = 1, 2, 

a(z*) = det[F[l](t, a;, ^),^[2](t, = 0. 

Therefore 



If n(X) > we have another analytic properties. Namely, a + (A) has an analytic continuation a(z) 
for z e C+ \ D + , and a - (A) has an analytic continuation o(z) for z G C_ \ D~ (a(z) = a*(z*)). 
The function b (A) has an analytic continuation b{z) for z 6 D~, and & + (A) has an analytic 
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continuation b(z) for z <G D + (and b*(z*) — b(z)). The spectral functions have the same asymptotics 
as above. 

If a(z) have zeroes Zj £ C+ \ D + (J = 1, 2, ....,£>) then 

afa) =det[Z[l](t,x,z j ),Y[2](t,x,z j )] =0. 
Hence vector-columns of the determinant are linear dependent: 



(3.5) Y[2}(t,x,z j )= lj Z[l}(t ) x,z j ), 7i = 



B( Zj ) _ A( Zj ) 



At the conjugated points z* £ C- \ D~ (j = 1, 2, 

a(z*) = det[F[l](^x,z;),Z[2](t ;a; ,z;)] = 0. 

Therefore 

(3.6) Y[l]{t,x,z*)= l3 ^[2]{t,x,z*l 7j = = = 7*- 

Being analytic in the domain D + the function b(z) can have zeroes at the points Zj € _D + (j = 
1, 2, g). Hence 



(3.7) Y[2](t,x,z 3 )=j j Z[2](t,x,z j ), ^ 



B( Zj ) A{z ) 



At the conjugated points z* € £)~ (j = 1,2, ....,q) 



(3.8) r[i](t,M;) =f i z[i](t,x,^), = = = 7 ;. 



4. Matrix Riemann-Hilbert problems 

In this section we give a reconstruction of the solution of the mixed problem to the MB equations 
in terms of the associated matrix Riemann-Hilbert problem. 
In the case n(A) < we put 

*(t,x,z)=\Z[l](t,x,z) Y ^^ Z) ^j ■ ^eC+, 

(4.1) 

= ( Y[l)(t,x,z) 
V a{z) 



Z[2]{t,x,z) , z£ 



Scattering relations (|3.ip gives: 

y-[l](t,a:,A) 
3 "(A) 

y + [2](t,a;,A) 
a+(A) 

where 



= Z-[l](t,x,X) -f~(X)Z-[2](t,x,X), 
= r + (X)Z+[l](t, x, A) + Z+[2](t, x, A), 



— (x) _ b (A) _ a-(A)i3(A)+/3+(A)^(A) + _ b+(A) _ a+(X)B(X) + (3 + (X)A(X) 

r{> 3"(A) 5 -(A)Z(A)-^-(A)S(A)' « + W «+(A)A(A) - P+{\)B(\) ' 

Using these relations we find that 

det $(t, x, X + iO) = det $(t, x, A - iO) = 1, Im A = 0. 
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Since $+(t, x, A) := $(t, x, X + iO) and a;, A) := $(£, x, A — iO) satisfy i-equation, we have 

<!_(£, x, A) = $+(t,x, A)J (x,A), 

where unimodular matrix 

W (r+wrn ^ Y+[2](0,x,X) Y 1 f Y-[l](p,x,\) 
J (x, A) = I Z+[l](0,x,A) 1 ^ a-(X) [2](0,x,A) 

All entries 

Z+[l](0,x,A) = W + (0,x,X)w + [l](0,X), 
F+[2](0,x,A) = VF+(0,x,A)$[2](0,A), 
F-[l](0,x,A) = T^-(0,x,A)$[l](0,A), 
Z-[2](0,x,X) = W-(0,x,\)w-[2](Q,\). 

of these matrices are known by initial and boundary conditions. Indeed, since 1^(0, x, A) 
w ± (x,X)(w ± {0,X))~ 1 then 



Z+[l](0,x,A) = tu+(x,A) 



01 ' 



F+[2](0,x,A) =w+{x,X) I &+(A) 



F-[l](0,x,A) = w-(x,X) 
Z-[2](0,x,A)=™-(x,A) (jj 



a+(A), 

a-jxy 

b (X)J ' 



The matrix Jo(x, A) can be written in the form: 

(4.2) J (x,A) = (S+(X))- 1 (w+(x,X)y 1 w~(x,X)S-(X) 

where 

6+(A) 



S + (X) = (l r+ [ X) y r+(A) 
S-W = (' % -(A)- ^ 



a+(A)' 



-r-(A) V' ^ v a "(A) 

are spectral matrices defined by initial £i n (t) and boundary £0(^)1 Po^A), AAo(a:, A) conditions. 
They are unimodular matrices. Matrix w^(x, A) is the Jost solution of x ± -equation (with t = 0): 

(4.3) -^-w ± (x,X) =(iXa 3 - iG ± (0,x,X) + H(0,x))w ± (x,X), 

ax 

u. ± (L,A)e- iL " ±(A)CT3 =1. 

It is defined by given boundary conditions only. The jump matrix Jo(x, A) can be rewritten in the 
form: 

J (x, A) = (K + (x, A))- 1 X-(x,A), 
where K (x, X) satisfies equation (I4.3[) under condition: 

K ± (L,X)=e iLri±as S ± (X). 
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The following reductions 

S~(X) =a 2 S+*{X)a 2 , 

W~ (x, A) =(72 w + * (x, A)(T 2 

K-{x,\) =a 2 K+*{x,X)a 2 

are fulfilled. These reductions are equivalent to the following ones: 

[^(A)]" 1 =[5"(A)]t, 
\w+{x,\)]-i=[w-{x,\)]\ 
[K+frX))- 1 =[K-(x,\)? 

where f means Hermitian conjugation. 

Lemma 4.1. For any x e K + and Ael the jump matrix 

J (x,X) = (K-(x,X)^K-(x,X) 

is positive defined. 

Proof. Indeed, scalar product 

(Jo(x,\)Z,Q = (K-(x,\)S,K-(x,\)Z) 

is positive for any (eC 2 ,^0. Suppose the contrary, i.e. there exists £ = €o{ x > X) such that 
K~(x, X)£o(x, A) = 0. Since detK~(x, A) = 1, vector function £o(x,X) is equal zero identically. 
Hence (J (x, A)£,£) > 0. □ 

Lemma 4.2. For any fixed x e R + the jump matrix Jq{x, X) has the following asymptotic behavior: 

J {x,X) = / + 0(A" 1 ), A^±oo. 
Proof. The matrix K := K~(x, A) satisfies equation: 

K' x = (iXa 3 + H(0, x) - iG_ (0, x, X))K 
Hermitian conjugation matrix satisfies the following one: 

Kl' = K\-iXa 3 + H\0,x) + iGl(0, x, A)) 
Since H^(0,x) — —H(0,x) and GL(0,a;A) = G+(0, x, A), the last equation takes the form: 

Kl' = Ki(-iXa 3 - fft( 0) x ) + iG+(0, x, A)). 
Take into account definition of Jo (a;, A) we find 

J' Qx = irf(G+(0,x,\) - G-(0,x,X))K = -^^-K^F(0, x, X)K. 

Since F(0,x, A) and if = if~(a:,A) arc bounded, and n(A) = o(A _1 ) as A — > ±oo, we have 
that J (x,A) = J(A) + o(A" 1 ) as A ^ ±oo. Then J(A) = J (0,A) + o(A- 1 ) and J o (0,A) = 
if+- 1 (0,A) J ft:-(0,A) =/ + 0(A- 1 ) as A^ ±oo. Hence J(A) = I + CKA- 1 ) as A -> ±oo. □ 
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For the problem on whole i-line the jump matrix Jq(x,X) can be found explicitly. Indeed, for 
t — > — oo we have 



a+(X) J \ a~(X) 
/ i\t-ix V+ (\) b+ ( x ) -ixt+ix V+ (x)\ ( e lA * +IX? '- (A) 



(4.4) 



V 



a+(A) 

p— iAt+ia;i7 + (A) 



/ 



b (A) 

iXt—ixr]- (A) iAi — i:c?7_ (A) 

V o-(A) 



W-V- PO) + y~ (\\ r + (X\ e ™(v+M-V-W) _ r +M\ e ix(7j+(A)— !7_(A))\ 



3 ia;(r/ + (A)-7j_(A)) 



where r (A) = r+ *(A). 
Now we put 

(4.5) 



M(t, x, z) = $(t, cc, z ) e i(«*-«iW)»-a ; z g C \ : 



where n(z) is defined by (|2 . 13|) . Then M(t, x, z) is a solution of the following Riemann-Hilbcrt 
problem RH tx : 

• M(t,x,z) is analytic (a{z) ^ 0) or meromorphic (a(z) has finite number of zeroes) in 
zeC\M and continuous up to the real X-axis; RHl 

• If a(zj) — a(z*) — 0,j = 1,2, ...,p then M(t,x,z) has poles at points z — Zj, z = z* 
(j = 1,2, ...,p), and corresponding residues satisfy relations: 



res z=Z] M[2](t,x,z) = m 3 e^ z ^ x ^ z ^ M[l]{t, x, Zj) 
res z=z . M[l](t, x, z) = m * e -^t-x n ( z *)) M[2](t, x, z*), 



RH2 
RHi 



where rrij = jj/a(zj), m* = "fj/a (z*), and numbers "/j, 7^ = 7* are defined in 
and [Jj\i: 

• M-{t,x,X) = M + (t,x,X)J{t,x,X), A G R, RHA 



(4.6) 



J(t,x,X) = e-' l{xt - x ^ +W)a3 J {x,X)e l{xt - x,1 - W) ' y3 - 



RH5 



• M(t, x,z) = 1 + 0(z _1 ), \z\ -> 00. 

Analytical properties follow from (14.11) . and residues relations arise from (|3.3[) - (|3.4[) . The jump 
matrix Jq(x, A) in (I4.6[) was described in (|4.2|) . (|4.3j) and Lemmas 4.1 and 4.2. The asymptotic 
behavior of M(t, x, z) follows from sections 2 and 3. 

For the problem on whole t-line the jump matrix degenerates to explicit one: 

J{t, X, A) =e -KA*-*U+(A))os j ( a . j A ) e i(At-x I7 _(A)) < r 3 



'l + |r+(A)| 2 e 2ia:( ' ? + (A: '~' ? - (A:)) _ r +(A)e~ 2iAt+2ia:J) +( A ) N 
-r+*(A)e 2iAt_2ia: ' ? -( A ) 1 
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7T1 

Here we have used (|4~4|) , (f372]) . Since rj± = ?y(A ± iO) = A - 7 (A) =p -J^PO where / (A) = 

^ n(s) (is 
p.v. then 

- oo s - A 4 

J(t,a:,A) = 7 + 0(e 7r "( A ^/ 2 ) -> 7, x -> +oo 

for a long attenuator (n(A) < 0), but it is unbounded (0(e 7rn ^^ x )) for a long amplifier (n(A) > 0). 
Therefore the jump matrix is exponentially closed to the identity matrix if n(X) < 0. Hence 
the main term of asymptotics contains solitons generated by discrete spectrum, while a contribu- 
tion of continuous spectrum is exponentially small. This justifies the phenomenon of self-induced 
transparency of attenuators. For the first time it was proved in [5]. 

In order to define a matrix RH problem for the mixed problem, we introduce another matrix 

Ht,x,z)=(z[l](t,x,z) m ^ ,Z) \ zeC + \D+, 

' Y[l](t,x,z) 

a(z) 
Y[2](t,x,z) 



Z[2](t,x,z)), z€C-\D~, 



b(z) 



Z[2](t,x,z) , zeD 



— r 1/ \ Y\i\(t,X,z)\ 

Z[l](t,x,z) 1 i> ' ' '- ), zeD-, 
b(z) J 

when n(A) > 0. By the same way, using residues relations (I3.5[l - (|3.8p . we obtain a meroniorphic 
matrix RH problem. But for the sake of simplicity we assume that a(z) ^ in the domain 
C+ \ D + and b(z) ^ in the domain D + , i.e. we consider a regular RH problem. In this case 
matrix M (t, x, z) — <!>(<, x, z )e 1 ( zt ~ ;CT '( z )) cr3 satisfies the next items. 

• M(t, x, z) is analytic z G C \ £, E := RU7 U7 and continuous up to the contour S; RRHl 

• M-(t,x,z) = M + (t,x,z)J(t,x,z), zeS, RRH2 

J(t, X, z) = e -m-*V+m)° S J ( X5 A J e i(At- jB i7_(A))a» ) A G K 

= e -i(»*-^W)"aj ^) e i(**-«>flW)»3 j z G 7 U 7, 

w/iere 

J (a;,z) =(^(a;,A))" 1 /f "(ic ) A), z = \eM. D =KnD_UD+, 



=(*-+(*, A))" 1 *^*, A), 



z = AeK\l 



■/<>(*) 



a(z) 

/ i 

\a(z) 

• det J(t, x, z) = 1 /or z G E; 

• M = Z+ 0(z _1 ), |z|^oo. 



a(z) 

1 J 

a(z)\ 
"&(z) 



/ 



z e 7, 



z g 7 



RRH3 
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a+(A) 



ffere iff (x, A) = w ± (x, X)Sq (A), 
S+(A) 



#±(x,A) = w ± (a;,A) 1 5* ± (A) and 

6 (A) 



5 "(A) 



S+(A) 



V6+(A) V 

a+(A) 
V° 1 / 



Vo 



i / 



s-(\) 



1 

6~ (A) 







V a- (A) / 



Note that matrix Jo(x, z) depends on x for z = A e I, and matrix Jq(z) is independent on x 
(and on t as well) for z € 7 C C+ and z£)C C_. Thus we prove 



Theorem 4.1. Let £(t,x), J\f(t,x,A) and p(t,x,X) be the solution of the mixed problem ljl.4\ )- 
Ji. 7[ ) to the Maxwell- Bloch equations \1. 1\) - FL3\) . Then there exists matrix M(t,x,z) which is the 
solution of the meromorphic Riemann-Hilbert problem (RHl)-(RHh) ifn(X) < or regular (under 
additional conditions a{z) ^ and b(z) ^ 0) RH problem (RRHl)-(RRHS) if n(X) > 0. Complex 
electric field envelope £(t,x) is defined by relation: 

(4.7) £ (t,x) = - lim 4izM 12 (t,x,z), 

and N{t, x, A) and p(t, x, A) can be found from linear equations by already known £(t, x). 

Proof. Formula (|4.7p follows from (|1.9|) and (RH5). Indeed, substituting (|4.5j) into equation 1(1.90 . 
we find 



Mt + iz[<T 3 ,M] + HM = 
m(t, x 



M(t,x,z) 



+ o(z~ 1 ), 



(4.8) 

Using (RH5) we put 
where 

m(t,x)= lim z(M(t,x, z) — I). 

z— )-oo 

This asymptotics and equation (|4.8p give 

H(t,x) = -i[a3,m(t,x)} 

and hence 

£(t,x) — — 4imi2 = — lim AizMi2(t,x, z). 



□ 



T7ms i/ie mixed problem on the finite interval < x < L (or half-line < x < oo) for the 
Maxwell-Bloch equations is linearizable completely. 

5. More general matrix Riemann-Hilbert problem 

Now we prove that any Riemann-Hilbert problem like RRH1-RRH3 generates a solution to the 
Maxwell-Bloch equations. From here and below we will consider more general construction. Let 
oriented contour £ contains real line K, circle T of sufficiently large radius and some finite arcs 
7j ■ U 7 j (j = 1, 2, ■■:,p), which are symmetric over the real line. Thus 

p 

s = Ruru [Jt,- utj-. 

3=1 
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If n(X) > then contour E includes additionally the closed oval 7 U 7 where Imrj(z) = 0. Such 
contours appear when we deal with periodic initial data or/and periodic boundary conditions. 
The large circle T allows to bypass difficulties connected with discrete spectrum and spectral 
singularities (sf.[21], [H]). Contour S has the following orientation: real line R is oriented from 
the left to the right, the circle T and oval 7U7 are oriented clock- wise, the arcs jj Ujj are oriented 
updown. Then the formulation of a regular matrix RH problem is as follows. 
Find 2x2 matrix M(t, x, z) such that 

• M(t, x, z) is analytic in z 6 C \ E and bounded up to the contour E; Rl 

• M_(t, x, z) = M+(i, x, z) J(t, x, z), z 6 E; R2 

• det J(t, x, z) = 1 for z G E; i?3 

• M(t,a;,z) = l + 0{z- x ), |z| 00. Ri 

Let contour E and jump matrix J(t, x, z) are satisfied the Schwartz reflection principle: 

• contour E is symmetric over real axis R, 

• J~ x (x, t, z) = J'{x, t, z*) for z £ E and Im z 7^ 0, 

where f and * are Hermitian and complex conjugations respectively. 
Furthermore, 

• jump matrix J(t,x, A) for A e R and iel has a positive definite real part and following 
asymptotic behavior: 

J(t,x,\) = I + 0(\- 1 ), A -> ±00. 

Theorem 5.1. Let jump matrix J{t,x,z) satisfies the Schwartz reflection principle, has a positive 
definite real part and I — J(t, x, .) G L 2 (T,) fl L°°(E). Then for any fixed t, x € R, the regular RH 
problem Rl, R2, i?3, i?4 /ias a unique solution M(t,x,z). 

Proof. Existence. Let t and a; be fixed. We look for the solution M(t, x, z) of the RH problem in 
the form 

(5.1) M(t,x,z)=I+± f P{t,x,s)[I-J{t, X ,s)\ _ 

27T1 J S — Z 

s 

The Cauchy integral (|5.ip provides all properties of the RH problem (cf.[2n]) if and only if the 
matrix Q{t, x, A) := P(t, x, X) — I satisfies the singular integral equation 

(5.2) Q(t,x,z)-K[Q](t,x,z)=R(t,x,z), z e E. 

The singular integral operator K, and the right hand side R(t, x, z) are as follows: 

fC[Q](t,x,z) := — / ds, 

27T1 J S — Z + 

s 

1 f I- J(t,x,s) , 
R(t,x,z) := — / v ' ds. 

2"7T1 J S — Z + 

s 

We consider this integral equation in the space £ 2 (E) 2x2 matrix complex valued functions 
Q(z) :— Q(t, x, z), z € S. The norm of Q € L 2 (T,) is given by 

1/2 / 2 \ 1/2 

\\Q\\m S ) = I / tr(Qt(z)Q(z))|dz| I = V / iq^W)! 3 !^! 
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The operator K. is defined by the jump matrix J(t, x, z) and the generalized function 

lim 



S — Z + z'^z,z'£side+ S — Z 

Furthermore, since the jump matrix J(i,x, A) has a positive definite real part, when A e I, then 
Theorem 9.3 from [21] (p. 984) guarantees the L 2 invertibility of the operator Id — K (Id is the 
identical operator). The function R(t, x, z) belongs to L 2 (Y>) because / — J(t,x,z) £ £ 2 (£) when 
z£S, and the Cauchy operator 

/(») _ /« . .. i l f(*l ds 



Zm J s — z + 2 27Ti J 



s — z 



is bounded in the space L 2 (Yj) |22| . Therefore, the singular integral equation (|5.2[) has a unique 
solution Q(t, x, z) £ L 2 (T,) for any fixed x, t € M., and formula (|5.ip gives the solution of the above 
RH problem. 

Uniqueness. A sketch of the proof is as follows. Since det J(t, x,z) = 1 one can find that 
det M(t,x, z) = 1 by repeating step by step the proof of the Theorem 7.18 from [20] (p. 194-198). 
Hence the matrix M~ l (t, x, z) exists, is analytic in z G C \ S and bounded up to the contour S. 
Let now suppose that there is another matrix M (t, x, z), which solves the given Riemann-Hilbert 
problem. Then 

M_ (i, x, z)MZ 1 (t, x, z) = M+(t, x, z)J(t, x, z) J _1 (i, x, z)M+ l (t, x, z) 

= M + (t,x,z)M+ l {t,x,z). 

since M(t,x,z) and M _1 (t, x, z) are bounded up to the contour S, end points and points of 
self-intersection are removable singularities. Hence the matrix M(t, x, z)M~ l (t, x, z) is analytic in 
z £ C and tends to identity matrix as z -> oo. By Liovilles's theorem M(t, x, z)M _1 (t, x, z) = I 
and therefore M(i, x, z) = M(t, x, z), i.e. the matrix M(t, x, z) is unique. □ 

Now we specialize the jump matrix J(t, x, z) in such a way that corresponding RH problem 
solution M(t, x, z) generates exactly the Maxwell-Bloch equations. Other specializations of the 
jump matrix J(t,x,z) give another nonlinear integrable equations. 

Theorem 5.2. Let inhomogeneous broadening n(X) be smooth and fast decreasing for AgR and 

oo oo 

rj(z) = z / — ^-ds, [ n(s)ds = ±1. 

4 J s — z J 

-co — oo 

Let z) := M(t, x, z)e~ 1 ( zt ~ x ' n ( z ^' J3 , where M(t, x, z) is the solution of the regular RH problem 

(Rl)-(RA) with the jump matrix J(t,x,z) such that 

' e -i(zt-x v (z))a 3 j( z } e i(zt-xri(z))<7 3j z £Y,\R (Imz^O), 

e -i(At-x^+(A))a 3 J( a . ) A ) e i(At-xn-(A))a 3) z = \ £ R; 

• J(z) is independent on t, x and satisfies the Schwartz reflection principle; 

• J(x, A) is independent on t and has a positive definite real part for any x £ M + and A £ M; 

• J(x, X) = I + 0(A- 1 ), A ±oo. 

If M(t, x, z) is absolutely continuous (smooth) in t andx, then <&(t,x,z) (z £ C±) satisfies AKNS 
equations: 

$t = -(iza 3 + H(t, x))$, <£ x = (izo-3 + H(t, x) - iG(t, x, z))$ 



J(t, x, z) 



20 



V. KOTLYAROV 



almost everywhere (point-wise) with respect to t and x. These equations become \1.9\) and \1.12)) 
as z = A ± iO. Matrix H(t, x) is given by 

H(t, x) = -i[<J3, m(t, x)], m(t,x) = — I (I + Q(t,x,z))(J(t,x,z) — I)dz, 

7T J 

where Q(t,x,z) is the unique solution of singular integral equation \5.2\) . Matrix F{t,x,X) is 

Hermitian and has the structure: F(t,x,X) — ( Jjh X ' ^!r}+^\\\- It * s reconstructed as 

\P \t,x,A) —AI{t,x,A)J 

the following jump: 

^^-F{t,x,X) = <$> x (t,x, z)^ 1 (t,x, z)| z =A+io - $ a; (t,a;,2;)$~ 1 (^x,z)| z= A-io- 
Proof. The matrix <£>(<, x, z) is analytic in z £ C \ £ and has the jump across S: 

$_(<, x, z) = $+(<, x, z)J(x, z), 
where J(x, z) is independent on t. This relation implies: 

^^(M,*) = d -^^^\t,x,z) 

for z £ E. This relation means that the matrix logarithmic derivative $t(i, x, z)<&~ l (t, x, z) is 
analytic (entire) in z £ C. Indeed, matrices M(t, x, z), M~ l (t, x, z) and derivative (in t) M t (t, x 7 z) 
are analytic in z £ C \ S and the Cauchy integral (|5.1I) gives the following asymptotic formulas: 

»r/, \ r , m±(t,x) _ 2 dM(t,x,z) dm ± (t,x)/dt _ 2 

M(t,x,z) = I-\ 1-0(2 ), ; = VO\z ), z £ C± 

z dt z 

as z — > oo. Hence 

$ t (t, x, z)^" 1 ^, x, z) = -iz(J3 + i[a 3 ,m + (t, x)] + 0(z _1 ) 

= — iz(73 + i[(T3, m_(t, a;)] + 0(z _1 ), z — > oo, 

where [A, B] := AB - BA and 

m-(t, x) = m+(t, x) = m(t, x) — — [ (I + Q{t 1 x, z))(J(t, x, z) — I)dz 

71" J 

a 

Since x, z)$~ l {t, x, z) has no jump across £ and M t (t,x,z), M _1 (i, x, z) are bounded up 
to the boundary, then the end points and points of self intersection of the contour £ are remov- 
able singularities for $f (t, x, z)$ _1 (i, x, z). Therefore, by Liouville's theorem, this derivative is a 
polynomial: 

U(z) := $ t (i, x, z)^" 1 ^, x, z) = -iza 3 - H(t, x), 
where H(t,x) := — i[cr 3 , m(t, x)] = ( ^ ^ ^q X ^\. Using the Schwartz symmetry reduction 

of the jump matrix J(t,x,z), we show fsf.|19]) that U(z) = 02U* (z*)<72, where 01 = q 

This reduction imply H(t,x) = —H'{x,£), i.e. q(t,x) = —p*(t,x) and we put q(t,x) := £(t,x)/2. 
Formula (15.11) gives an integral representation for the electric field envelope 

(5.3) £(t,x) = ^-J {[I + Q(t,x,z)][J(t,x,z) - I}) n dz. 

s 

through the solution Q(t,x,z) of the singular integral equation (|5.2j) . Thus $(t, x, z) satisfies 
equation (|1.9|) . and a scalar function £(t,x) is defined by (|5.3|) . 
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In contrast with previous case matrix logarithmic derivative $ x (t, x, (t, x, z) is analytic in 

z G C± only. Indeed, since matrix J(z) is independent on t and x for z £ E\R, then this logarithmic 
derivative is continuous across the contour E \ R, while it is not continuous across the real line 
because J(x, A) depends on x. End points of the contour E are removable singularities because 
the matrices M x (t 7 x, z) and M _1 (i, x, z) are bounded up the conjugation contour E. Further, the 
asymptotic behavior at infinity gives 

= izer 3 + H (t, x) + 0(z 1 ), z e C±, z -> oo. 

Therefore we find that & x (t, x, 2;)$~ 1 (t, x, z) — iz<r 3 — H(t, x) is represented as the Cauchy integral: 

oo 

1 f F(t x s)n(s) 

^ x {t,x,z)^~ l {t,x,z) - iz<7 3 - H(t,x) = — / v ' ' ; y ' ds, z R, 

4i J s — z 

— oo 

where F(t, x, A) is some matrix, and factor l/4i is chosen for a convenience. Due to reflection 
(symmetry) property of the jump matrix we find that F(t, x, A) is Hermitian. Since tr($ x (t, x, A ± 
iO)$ _1 (i,a;,A±iO)) = (det x, A ± iO))!, = and tra 3 = trH(t,x) = then trF(t, x, A) = and, 
hence, F(t,x,X) has the structure: 



F(t,x,A):=f^J*' a! '^ P A ( A *f' A 1,V 

v 7 V/ 9 (*' X ' A ) -M{t,x,\)) 



Thus we see that the matrix $(t, x, z) satisfies two differential equations: 

(5.4) $ t = ?7(i, x, z)®, U(t,x,z) = -iza 3 -H(t,x) 

(5.5) $ x = V(t,x,z)&, V(t,x,z)=iza3 + H(t,x) — iG(t,x,z), 
where 



G(t,x,z) = - / ds, zf 

4 J s - z 



For real z = A e R we have two differential in x equations: 

(5.6) § x = V±(i,x,A)$, V±(f,x,A) = iAcr 3 + H(t, x) - iG±(t, x, A), 

where G±(i, x, A) := G(t, x, A ± iO). In particular, the last equations give a reconstruction of the 
function F(t, x, A): 

m W p ^ ^ ^ = x ^ x + io)^- 1 ^, x, A + iO) - $ x (t, x, A - iO)$ _1 (t, x, A - iO), 

where $(t, x, z) = M(t, x, z)e~^ zt ~ xr, ^ c ' 3 , and M{t,x,z) is the solution of the regular RH problem 
(iil)-(-R4). The compatibility condition (& x t{t, x, A ± iO) = &tx{t, x, A ± iO)) gives the identity in 
A: 

U x (t, x, A) - V±{t, x, A) + [U(t, x, A), V ± (t, x, A)] = 0. 
This identity is equivalent to 



H t (t,x) +H x (t,x) - i J [cr 3 , F(t,x, s)}n(s)ds = 

— oo 

o 

• (F t (t, a;, s) + [isCT 3 + H(t, x), F(t, x, s) 



oo 

i 



s - A T iO 



n(s)ds 
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and it is possible if and only if the left and right hand sides are equal zero, i.e. 

oo 

H t (t,x) + H x (t,x) - - J [a 3 ,F(t,x,s)]n(s)ds =0 

— OO 

F t (t,x,X) + [iXa 3 + H(t,x),F(t,x,X)] = 0. 

These matrix equations are equivalent to the MB equations (|1.1[) - (|1.3|) . Thus we proved that 
the matrices $>(t,x, X ± iO) satisfy equations (15.41) and (15.61) . which coincide with AKNS system 
(|1.9p and (| 1 . 12|) . and scalar functions £(t,x), J\f(t,x,X) and p(t,x,X) satisfy the Maxwell-Bloch 
equations (fTT ]) -(fI3 1) . □ 



6. Conclusions 

It is proved that the mixed problem (11.4[) - (|1.7p to the Maxwell-Bloch equations (jl.l|| - (|1.3|l 
is linearizable completely by using the matrix Riemann-Hilbert problem RH1 - RH5 or RRH1 - 
RRH3. More general Riemann-Hilbert problem Rl - R3 generates a solution to the Maxwell-Bloch 
equations if conjugation contour and jump matrix are satisfied the Schwartz reflection principal 
and some special restrictions. Among generated solutions there are: solutions defined for t £ K and 
x £ M+ and studied in [5], [6], solutions to the mixed problem ([LT jl -lfO ]) . ([L3 ]) -(fI7T |) (t,x G K+) 
with decreasing or periodic input pulse £{t 1 0) and different initial functions £(0, x), Af(0, x), p(0, x), 
etc. The kind of solutions is defined by the specialization of the conjugation contour and the jump 
matrix. Suggested matrix RH problems will be useful for studying the long time/long distance 
asymptotic behavior of solutions to the MB equations using the Deift-Zhou method of steepest 
decent. 
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